Subthermal linewidths in photoassociation spectra of cold alkaline earth atoms 
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Narrow s-wave features with subthermal widths are predicted for the 1 Yl g photoassociation spectra 
of cold alkaline earth atoms. The phenomenon is explained by numerical and analytical calculations. 
These show that only a small subthermal range of collision energies near threshold contributes to 
the s-wave features that are excited when the atoms are very far apart. The resonances survive 
thermal averaging, and may be detectable for Ca cooled near the Doppler cooling temperature of 
the 4 1 P^4 1 S laser cooling transition. 



O 

Oh! 



in 

o 

o 



CZ3 

O 

Oh- 



X 



I. INTRODUCTION 

Photoassociation spectroscopy has become a very pow- 
erful tool for studying the collision physics of laser cooled 
and trapped atoms m . The conventional wisdom is that 
the linewidth of individual molecular levels in the pho- 
toassociation spectra of laser cooled atoms is due to the 
natural linewidth plus thermal broadening on the order 
fcflT, where ks is the Boltzmann constant and T is the 
temperature. Thus, we would not expect photoassocia- 
tion lines to be much smaller than fc^T in width [§-0. 
However, we demonstrate subthermal linewidths for a 
special case of photoassociation at very long range to 
an excited vibrational level v with a small natural decay 
width T v . In this case the s-wave vibrational features are 
very narrow at low T (e.g., < 10 //K), where fcgT -C T v . 
Surprisingly, such features can remain narrow even at 
much higher T (e.g., ss 1 mK), where fcgT 3> r„. These 
subthermal linewidths are a consequence of having only 
a narrow range of low collision energies that contribute 
to the thermally averaged photoassociation spectrum. 

Narrow features are possible in 1 H g photoassociation 
trap loss spectra at small detuning of alkaline earth atoms 
in a magneto-optic trap (MOT) j|. For Ca there is a 
good chance that the sharp s-wave features will stand 
out on a background comprised of broader peaks from 
the higher partial waves of the 1 H g spectrum and the 
broad features of the £« spectrum, even around the 
Doppler cooling temperature Tn (To = 0.83 mK for 
the Ca 4 1 P^4 1 S cooling transition). Photoassociation 
spectroscopy in a Ca MOT has been reported || , but in 
that experiment the photoassociating laser was detuned 
far from atomic resonance (about 780 atomic linewidths 
r at ), whereas our features are predicted for small detun- 
ings (< 25 T at ). 



II. THEORY OF SUBTHERMAL LINE SHAPES 
A. Trap loss spectrum 

In Ref. [pi we outlined the numerical and analytical 
models used here. The numerical photoassociation trap 
loss spectrum is obtained from a fully quantum mechan- 
ical three-channel model, where the time-independent 
Schrodinger equation is solved with a complex potential 
to represent spontaneous decay from the excited state. 
We repeat only the essential analytical formulas here, 
and concentrate on the explanation of the phenomenon. 
For a more detailed description of different aspects of 
cold photoassociation collisions of alkaline earth atoms, 
see Ref. §. 

Figure ID shows a high resolution scan of our calculated 
state change (SC) trap loss spectra of cold Ca atoms at 
three different temperatures for a single vibrational level 



of the 1 II 9 state. 



There are a number of 1 Tl g isolated vi- 



brational resonance features in the range of detuning A 
from 1 to 25 T at to the red of atomic resonance. Our ex- 
ample is typical of these features. Note that the features 
labeled A and B in the figure are both narrower than 
ksTo- Although the transition between the 1 T, g ground 
state and the 1 n ff excited state is forbidden at short and 
intermediate internuclear distances R, it becomes allowed 
at long range due to relativistic retardation effects P,10[ . 
The trap loss collision of the two cold atoms proceeds 
via excitation at a long range Condon point Re to the 
1 Tl g state (the difference of ground and excited molecu- 
lar potential energy curves equals the photon energy at 
Re)- Once in the 1 Ii g state the atoms are accelerated 
towards short range. The survival probability in moving 
from long to short range on the ^g state is close to unity 
due to the small decay rate (T(R) — ► for R — ► 0). At 
short range a SC may occur due to spin-orbit coupling 



to a lower lying state correlating to atomic 1 S + -"-D, 3 D 
or 3 P states. After SC to these channels, modeled here 
by a single effective channel, the atoms will be lost from 
the trap due to the large gain in kinetic energy. 

The photoassociation spectrum in Fig. tt\ is the ther- 
mally averaged loss rate coefficient flf]] : 

k(t,a)= v ( M+ i^/MiMlE 
i evsn ,j v \ k 

= £ (K(e,A,£J)), (f) 

ieven,J 

Here e = H 2 k 2 /(2fj,) is the collision energy at momentum 
hk for reduced mass fi, £ is the ground state partial wave 
quantum number (0, 2,4,... for identical Group II spin- 
less bosons), J = £,£ ± 1 is the excited state rotational 
quantum number, and A is the detuning from the atomic 
resonance of the photoassociating laser. S pg (e, A,£J) 
is the S'-matrix element for the transition between the 
ground state g and the SC channel p via the excited state 
e. The brackets (• • •) imply a thermal average over a 
Maxwellian energy distribution: 



x l ' 2 e 



x {- ■ -)dx 



(2) 



where x — e/iksT). As we show below, a consequence 
of the excitation at a very large Rc is that only a small 
range of energies, much less than fc^T when T is near Td, 
contributes to the thermal average integral in Eq. (0), 
especially for s-waves. Consequently, averaging does not 
introduce much additional broadening. 



B. Analytic theory 

An analytic interpretation can be given for the ori- 
gin of the subthermal linewidths. When the spacing hv v 
between vibrational levels v is much larger than their 
total width r„, i.e., the vibrational resonances are non- 
overlapping, then \S pg \ 2 is given by an isolated Breit- 
Wigner resonance scattering formula for photoassocia- 
tion lines @,|,|: 



\S pg (e,A,£J)\ 2 



*■ vp^ vg\£-) %-J ) 



[e-£ v (A,J)} 2 + (T v /2) 2 ' 



(3) 



The total width r„ is the sum of the decay widths into 
the SC (T vp ) and the ground state (r„ fl ) channels and 
the radiative decay rate (r„, r ad); and £ V (A, J) = A — 
[s v (J) + s v (J)] is the detuning-dependent position of the 
vibrational level v J in the molecule-field picture relative 
to the ground state separated atom energy. The level 
shift s v (J) due to the laser- induced coupling Q is small 
for our case. When A = e v (J) + s v (J), then £ V (A, J) = 
and the vibrational level is in exact resonance with 
colliding atoms with zero kinetic energy. 
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FIG. 1. Single vibrational feature of the Ca 1 Yl g photoas- 
sociation spectrum resulting in SC trap loss. The A/T at value 
where the peak appears depends on the model potentials in 
Ref. pi. It may appear at a different A/F a t in an actual ex- 
periment. The smooth background of about 2 units due to the 
1 T, g — > 1 E„ transition || is not shown. The Doppler cooling 
temperature for cooling on the 4 X P atomic state is 830 fiK, 
and the recoil temperature is 2.7 /iK. The low temperature 
limit for the feature labeled A is a Lorentzian line centered at 
the resonance position A/T a t = 11.37 with a natural width of 
1.2 MHz. The peak of the A feature, which clearly has a sub- 
thermal width for T = 830 /iK, is normalized to unity for the 
three different temperatures. The calculated peak trap loss 
rate coefficients for a 1 mW/cm 2 laser are 2 x 10~ 13 cm 3 /s 
(830 ,uK), 2 x 10~ 12 cm 3 /s (83 /jK), and 8 x 10 -12 cm 3 /s (8.3 
/iK), respectively. 



In the reflection approximation T vg {e,£J) is propor- 
tional to the square of the ground state wavefunction <p g 
at the Condon point (R c ) iP@: 



-L vg \&i %-J ) 



2Trh Vv \V eg (R c ,U)\< 
D c 



\cP g (eJ,R c )\ 2 . (4) 



Here V eg (Rc,iJ)) is the laser- induced coupling, Dc is 
the slope difference of the ground and excited state po- 
tentials at Rc, and v v is the vibrational frequency for 
level v. \V eg \ 2 is linear in laser intensity / for our as- 
sumed weak-field case. Approximating the ground state 
wavefunction by its low-energy asymptotic form gives for 
s-waves 



\ct> g (e,Q,Rc)\' 



2/j. sin 2 k(R c - A ) 
rf? k : 



(5) 



where Aq is the scattering length of the ground state 
potential. For higher partial waves (£ > 0) Ej] , 



\Me,£,Rc)? 



2fi z 2 c \ 3l (z c )\< 



■Kff 



k 



(6) 



where ji(zc) is the spherical Bessel function, and zc — 
kRc- The normal scattering phase shift does not appear 



in Eq. (||) since it is vanishingly small near threshold for 
the higher partial waves JJJ. We define here the near 
threshold range of collision energies to be the range for 
which Eqs. (|5j) and (|6|) are good approximations. 
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FIG. 2. Thermally averaged Ca SC trap loss rate for T = 
0.83 mK and I — 1 mW/cm 2 from the analytic line shape 
model, Eqs. (pi), (|5|) and (pi), compared to the full quantum 
numerical result, including contributions from £ < 6. The 
separate contributions are shown for the $-, d-, and g-waves, 
with the dominant contribution to each peak labeled by U. 

Figure shows the thermally averaged Ca trap loss 
spectrum K(T, A) at To = 0.83 mK obtained by insert- 
ing Eq. (g) into Eq. (g) using Eq. (g) or (|). The analytic 
model takes a few input parameters from the numerical 
model @] : the scattering length for the ground state po- 
tential (the actual value is unknown, but the model value 
is 67 a ), T vp /h w 0.3 MHz < T v>rad /h = 0.8 MHz, and 
the position e v (J)+s v (J) of the vibrational levels for each 
J. The good agreement with the quantum numerical cal- 
culations indicates the quality of the analytic model. 

Figure |2| also shows the individual contributions from 
the s-, d-, and g- waves (£ — 0,2,4). The analytic for- 
mulas also compare very well with the details of these 
individual features in the numerical calculation (compar- 
ison not shown) . We will concentrate on the overall sub- 
thermal features labeled A and B. The A feature is a 
sharp si line sitting on a background due to d and g 
lines, whereas the B feature takes its relative sharpness 
from a d3 line sitting on the d2 and g4 background. The 
dl feature contributes the shoulder on the left of the A 
peak. 



C. Origin of subthermal features 

Figure g shows K (e, A, £J) obtained from the analytic 
model for the £J = si, d2, and d3 features in Fig. @. 
The slanted maxima follow the lines of exact resonance 
where e = £ V (A,J). Figure || also shows a cut of the si 
and d2 K(e,A,£J) at the fixed detuning A/T at = 11.37 



where £ v (A, J) = for the si resonance (d3 has a similar 
energy variation as d2, centered at a different detuning). 
The variation of the integrand K(e,A,£J) in Eq. (0) as 
a function of e and A provides an explanation for the 
subthermal features. Since K(e, A,£J) peaks in a small 
range of e <C kriTo , only a small range of collision en- 
ergies, e/ks ~ 0.1 mK for the s-wave and 0.2 mK for 
the d-wave, contributes to the width of the feature (for 
comparison, T Vira d/kB — 0.04 mK). Consequently, the 
s-wave peak A broadens only slightly when the temper- 
ature increases from T = 0.083 mK to the Doppler limit 
To = 0.83 mK, as seen in Fig. [j]. Even the d-wave fea- 
ture B remains subthermal, although broader than the 
s-wave feature. 
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FIG. 3. The lower panel shows the analytic K(e, A,^J) in 
atomic units (1 a.u. = 6.126 x 10 -9 cm 3 /s) versus e and 
A for U = si, dl, and d3 (/ = 1 mW/cm 2 ). The vertical 
and horizontal axes are set to have the same energy scale in 
common units (ks/h — 20.8 MHz/mK) so that the dash-dot 
lines of exact resonance, e = £ V (A,J), have a 45 degree slant. 
The dotted contour lines start at 0.000025 and increase in 
steps of 0.000025. The solid contour lines start at 0.0002 
(bold) and increase in steps of 0.0002. The upper panel shows 
a cut of K(e, A, £J) for £J — si, d2 along the dashed vertical 
line of constant detuning indicated on the lower panel. The 
analytic threshold range e t h is indicated. 

The variation of K{e,A,£J) is a consequence of the 
near-threshold resonance form of S pg (e,A,£J), Eq. (||). 
The resonant denominator causes the largest contribu- 
tion to K (e, A, £J) to come from energies near the slanted 
line of exact resonance, e = £ V (A,J), indicated in the 
lower panel of Fig. 0. On the other hand, the term 



T vg (e,£J) [Eq. (§)], proportional to |</> g (M, Rc)\ 2 /k, is 
strongly influenced by the near-threshold properties of 
<p g (e,£,Rc)- We may distinguish two regimes: k <C kth 
where \<j> g (e,£, Rc)\ 2 /k ex k 21 , and k » /c t /j, where 
|^> 9 (e,^, i?c)| 2 /^ oscillates with an amplitude decreasing 
as 1/k 2 oc 1/e. Thus, the integrand K(e, A,£J) for s- 
waves approaches a constant value for k <C fcth and drops 
off rapidly and oscillates when k 3> kth ■ This variation is 
evident in the upper panel of Fig. || Using Eqs. (|J) and 
@, we estimate fc t h, from k t h\Rc — A \ = 7r /2 for s-waves 
and k t hRc — Z\{€) for £ > 0, where the first maximum 
in j'^(z) for positive argument is at z = zi(^). Taking i?c 
= 513 ao and the arbitrary model value v4_o = 67 ao for 
our case gives Sth/ks — [frkth) 2 /(2^ifes) = 0.05 mK for 
the s-wave and 0.18 mK for the d-wave. Thus, the large 
Re leads to the small value for kth and Sth, and conse- 
quently to the subthermal linewidth. For energies higher 
than eth, K(e 7 A, £J) has a node at e where <p g (e 7 £, Re) 
has a node at Re, e.g., at e/kg — 0.21 mK for the s-wave 
and 0.55 mK for the d-wave. A second maximum in the 
ground state wavefunction appears at higher energy, 0.48 
mK in the case of the s-wave. This is the origin of the 
maximum in the red wing of the si feature labeled si' in 
Figs. | and |. 

Figure [| illustrates the dependence of the si feature 
on the unknown ground state scattering length Ao- Since 
hh = tt/(2|-Rc — A) I), hh decreases for A < 0, resulting 
in narrower thermally averaged lines. However, if Aq is 
positive and near Re so that \Re — Aq\ becomes small, 
then k t h increases and the narrow peak broadens and 
flattens, so that it may no longer stand out. 



III. CONCLUSION 
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FIG. 4. Variation of the thermally averaged si feature at 
0.83 mK for three different ground state scattering lengths 
(Ao — 67 ao for Figs. yj-g). Thermal averaging causes the 
peak to shift from the exact resonance position A/P a t = 11.37 
(where £ v — 0) depending on the value of Aq. The si fea- 
ture broadens and shifts considerably when Ao becomes close 
to Re- In this case, e t h increases and a much larger range 
of collision energies contribute to the thermal average. The 
feature has a linewidth closer to the natural linewidth of 1.2 
MHz and shows more prominent subsidiary maxima when the 
scattering length is negative. In this case, eth decreases and 
a smaller range of energies contribute to the thermal average. 



small, < k <C |.Ao| _1 , it would be interesting to see if 
subthermal lineshapes are possible for very narrow res- 
onances with widths <C fcsT when \A \ is large enough 
that h 2 l(2[iAl) < k B T. 



We predict that subthermal line shapes should appear 
in high resolution photoassociation spectra of the ^g 
state of Ca dimer near the 1 P^- *S laser cooling transi- 
tion. Such features will be hard to see for Mg, where weak 
1 n s lines are obscured by a large 1 S M background M. 
Subthermal lines may be less prominent for Sr or Ba 
because of additional predissociation broadening and a 
higher density of states that combine to give broader and 
more blended ^g features. 

Subthermal linewidth of scattering resonances are pos- 
sible when the contributions to the line shape from the 
relevant S-matrix elements is restricted to very low col- 
lision energies below the range of thermal energies ksT. 
In our present study, this is a consequence of the very 
large Condon points associated with the transitions. It 
would be useful to extend this analysis to other pho- 
toassociativc transitions or magnetically-induced Fesh- 
bach resonances Eq]. This would be most interesting in 
the case of a large s-wave scattering length, that is, when 
\Aq\ 3> xo, where £ = ^(^^Cq/H 2 ) 1 ^ 4 is a characteristic 
length scale for a van der Waals potential with disper- 
sion coefficient Cq M. Since the near threshold range is 
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